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Abstract 

> ■ 

^ ■ We compare the Wronskian method (WM) and the Schrodinger eigenvalue march 

or canonical function method (SEM-CFM) for the calculation of the energies and 
eigenfunctions of the Schrodinger equation. The Wronskians between linearly inde- 

> , 

' pendent solutions of the Schrodinger equation provide a rigorous basis for some of 

' the assumptions of the SEM-CFM, like, for example, the concept of "saturation". 

vn 

We compare the performance of both approaches on a simple one-dimensional model 
and suggest that taking into account the asymptotic behaviour of the wavefunction 
(as is already done in the WM) may make the SEM-CFM more efficient from a 
' numerical point of view. 



1 Introduction 



In a recent paper Tannous and Langlois [T] proposed the Schrodinger eigen- 
value marcli (SEM) for the calculation of eigenvalues of the Schrodinger equa- 
tion. The SEM was developed by Kobeissi [2] under the name of canonical 
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function method (CFM) and even Tannous and Langlois [3] and Tannous et 
al used to call it that way before renaming it. According to the authors, 
the SEM compares favourably to the well-known Numerov and shooting meth- 
ods [1] . Leubner [5j argued that the underlying idea of the CFM is very old and 
can be found in standard textbooks on differential equations as the reduction 
of two-point boundary value problems to initial value problems. In addition to 
it, he also proved the equivalence between the CFM and the widely used shoot- 
ing method. In a rather late reply, Kobeissi |6] tried to show that the CFM is 
considerably more accurate that the shooting method and other algorithms. 
However, Tellinghuisen confirmed Leubner's conclusions and showed that 
the apparent numerical advantage found by Kobeissi [6] was merely due to dif- 
ferent numerical precision in the results compared by this author. Leubner's 
and Tellinghuisen's papers have been utterly omitted in all later applications 
of the CFM [4J. It is worth quoting Tellinghuisen's comment that 'Although 
the CFM is fundamentally less efficient than the Cooley algorithm, it does offer 
some advantages in practical applications, as it avoids unnecessary integration 
in the nonclassical region'. 

We have recently discussed a method for the calculation of bound states and 
transmission probabilities for one-dimensional wells and barriers p-[TU]. It is 
based on the fact that the coefficients of the linearly independent asymp- 
totic contributions to the wavefunction can be easily expressed in terms of 
Wronskians. Those papers were mainly focused on the pedagogical value of 
the Wronskian method (WM) that is already known since long ago (see the 
references cited there P4T0]). Since the WM is also based on the so called 
canonical functions [2]-|ll|6] (or normalized solutions pLj) we think that it may 
be interesting and fruitful from a pedagogical point of view to compare it with 
the SEM (or CFM). In Sec. [2] we outline the main ideas behind the WM. In 
Sec. [3] we discuss the SEM-CFM and provide a rigorous basis for some of 
its equations by means of the WM. In Sec. IHwe test, discuss and verify the 



2 



general results by means of simple examples. Finally, in Sec. [5] we summarize 
the main results and draw conclusions. 



2 The Wronskian method 

In order to introduce the main ideas of the WM we consider the second-order 
differential equation 

L{y) = y"{x) + Q{x)y{x) = Q (1) 

If yi and y2 are two solutions to this equation then we have 

yiL{y2) - y2L{y,) = ^W{y,, y^) = (2) 

where 

W{yi,y2) = yiy'2-y2y'i (3) 

is the Wronskian (or Wronskian determinant). If yi and y2 are linearly inde- 
pendent then the Wronskian ([3]) is a nonzero constant. For practical purposes 
it is convenient to choose two solutions C{x) and S{x) that satisfy 

C{xo) = S\x,) = 1, C"(xo) = ^(xo) = (4) 
at a given point Xq, so that W{C, S) = 1 for all x. 

In this paper we are interested in the dimensionless Schrodinger equation 

- + v{x)^{x) = e^{x) (5) 

where e and v{x) are the dimensionless energy and potential, respectively. In 
earlier papers we have discussed a systematic way of converting the standard 
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Schrodinger equation into its dimensionless form [814T0] and we do not repeat 
it here. In order to facilitate the discussion below we assume that — oo < x < 
oo, keeping in mind that the main results can be extended to other cases if 
necessary. If C{x) and S{x) are two solutions satisfying (|4]) then we can write 
a general solution as 

ifix) = A2C{x) + B2S{x) (6) 

where, obviously, 

A2 = ^{x^) = W{v,S) 

B2 = ^\xo) = W{C,^) (7) 

For notation simplicity we write C{x) and S{x) instead of the more detailed 
expressions C(e, Xo,a;) and S'(e, Xo,x), respectively, that explicitly indicate the 
dependence of the linearly independent solutions on the dimensionless energy 
e and the chosen coordinate point xq. 

It is well known that for an arbitrary value of the dimensionless energy e the 
wavefunction behaves asymptotically as 



AiLc{x) + BiL(i{x), X —7- — oo 
A^Rcix) + BsRdix), X ^ oo 



where L and R stand for left and right and c and d for convergent and diver- 
gent, respectively. It means that, for arbitrary e, the wavefunction is a linear 
combination of a convergent and a divergent function when |x| — i- oo. If, for a 
particular value of e, Bi = B^ = then the resulting wavefunction is square 
integrable. This condition determines the energies of the discrete spectrum. 

It follows from the well known properties of the Wronskians p-[TO] (and ref- 
erences therein) that 
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B,W{R,, R,)+ = A^WiR,, C)+ + B^WiR,, S)+ (9) 

where the subscripts — and + indicate that the Wronskians are calculated in 
the limits x — ?■ — oo and x — )■ oo, respectively (where they become constants). 
Therefore, when i?i = i?3 = we have a linear homogeneous system of two 
equations with two unknowns: A2 and B2. There will be nontrivial solutions 
provided that its determinant vanishes 

W{L,, C).W{R,, S)+ - W{R,, C)+W{L,, S)^ = (10) 

The roots of this equation e„, n = 0,1, . . ., are the energies of the bound states. 
When the potential is parity invariant 

v{—x) = v{x) (11) 

and a;o = then C{x) and S{x) are even and odd functions, respectively. In 
this case we have 

W{L„S). = W{R„ S)+ 

W{L„C)^ = ~WiR„C)+ (12) 

and the determinant (fTUjl takes a simpler form: W{Rc, C)+W{Rc, S)^ = 0. 
We appreciate that the even and odd solutions are clearly separate and their 
eigenvalues are given by 

W{R„C)+ = 

W{R„S)^ = (13) 
respectively. Besides, we need to consider only the interval < a; < 00. 
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3 The canonical— function method or Schrodinger eigenvalue march 



In order to compare the results of Sec. [3] with the SEM-CFM [lr4j we sim- 
ply note the following equivalence between the main functions ip{x) — )■ y{x), 
C{x) — 7- a{x) and S{x) — )■ f3{x). It follows from equation ([7j) that 

ip{xo) W{ip,S) ipS'-if'S ^ ' 

In the SEM-CFM one defines 

and obtains the eigenvalues from the roots of P-H] 

F(e) = /+(e)-/_(e) = (16) 

Note that Eqs. (fT5l) are identical to Eqs (4) of Tannous and Langlois [1]. 
Therefore, the WM is quite similar to SEM except for some slight differences 
that we will discuss later on. 

For concreteness we assume that the acceptable solutions to the dimensionless 
Schrodinger equation ([5]) satisfy the boundary conditions 

lim ^[x) = (17) 

|a;|^oo 

In such a case it is customary to simplify the SEM-CFM equations (fT5|) as 
follows [M] 



Ue)= lim ^^''^ 



x—^—oo 



Six) 



hie) = lim (18) 
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In order to provide a rigorous proof for these equations we take into account 
that 



C{x) = aiLc(x) + biLd{x), x — )■ — oo 
S{x) = a[Lc{x) + b[L(i{x), x -> — oo 
C{x) = a^Rc{x) + h^Rd{x), x — )■ oo 

S{x) = a'^Rc{x) + b'r^Rdix), X oo (19) 
so that 

C{x) hi 



l_ = hm 



S{x) b[ 



1 r ^3 .... 

We can thus give a precise meaning to the word "saturation" often used in 
connection with the SEM-CFM [Il[3lll]. It simply points to the obvious fact 
that the divergent functions dominate when |x| is sufficiently large and l± tend 
to the ratios of their coefficients in the expansion of the functions C{x) and 
S{x). Besides, the WM gives us expressions for the coefficients in equation 
(fT9i) and their ratios read 



h_W{C,L,)_ 

(21) 



b[ W{S,L,)_ 
h W{C,R,). 



Consequently equation (fT6|l leads to (fTOj) that has been rigorously proved in 
Sec. [21 Obviously, the right-hand sides of equation ( 12U]) do not change if we 
substitute C'{x)/S'{x) for C{x)/S{x) as pointed out by Kobeissi [2] without 
giving a rigorous proof. More precisely, if we substitute any pair of linearly 
independent solutions for C{x) and S{x) in equations f l20|) and f l2T|) we should 
obtain a similar result. It does not mean that all the pairs of solutions will 
be equally efficient from a numerical point of view. We simply want to point 
out that the WM provides a rigorous proof for the main equations commonly 
used in the SEM-CFM. 
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If we require that (p{xl) = v{xr) = for ^ Xq <^ x/j then we have an 
homogeneous system of two equations with two unknowns 



A2C{xl) + B2S{xl) = Q 
A2C{xr) + B2S{xr) = Q 



(22) 



that has nontrivial solutions only if 



C{xL)S{xn) - C{xr)S{xl) = 



(23) 



This determinantal equation, derived earlier by Leubner [S], is equivalent to 



Z+(e) — /-(e) = since /+(e) = C{xr)/S{xr) and /-(e) = C{xl) / S{xl) in the 



right and left asymptotic regions, respectively. 

In the case of a symmetric potential and Xq = we have C{—x) = C{x) and 
S{—x) = —S{x) and the approximate eigenvalues are roots of the simpler 
equation C{xr)S{xr) = 0; that is to say C{xr) = or S{xr) = for the even 
or odd solutions, respectively. 

4 Examples 

As a first illustrative example, Tannous and Langlois [T] considered a particle 
of mass m in a box with impenetrable walls at x = and x = a 



as a model for an electron in a one-dimensional metallic rod of finite length. 
The electron moves freely inside {V{x) = 0) and cannot escape from the rod 
{V{x) = oo if X < or X > a). They applied the SEM and obtained the 
well known energies. However, their choice of the model parameter a is rather 
atypical. Here, on the other hand, we transform the Schrodinger equation f l24|) 
into a dimensionless eigenvalue equation by means of the change of variables 



2m 



riX) = Eij{X), ^(0)=^(a)=0 



(24) 
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X = ax and ^{x) = ^/a■^p{ax) that leads to (f"{x) = —2e(f{x), where e = 
ma^E/h'^. Note that the boundary conditions for the dimensionless solutions 
are <^{0) = <^{1) = 0. 

Tannous and Langlois mention the problem of defining proper self-adjoint ex- 
tensions of the operators for the infinite-well potential [Hj. This mathematical 
subtlety is of great importance in the discussion of physical observables but it 
is not an issue with regard to the calculation of the eigenvalues by means of a 
numerical method like the SEM-CFM. 

The particle in an infinite square well is suitable for the pedagogical analysis 
of the performance of the shooting methods. If we consider a set of discrete 
coordinate points Xj = jh, j = 0,1, N such that x^ = Nh = 1 and define 
the approximate finite-difference first and second derivatives 

^{x + h)- ipjx - h) h^^"'{x) 
6h^{x) = — = ^{x) + + . . . 

ipjx + 2h) - 2^{x) + ifjx - 2h) h^p^^ 
dhVix) = ^ = ip {x) + + . . . (25) 

then the Schrodinger equation becomes a three-term recurrence relation 



+2 + {Shh - 2) cpj + ipj_2 = (26) 



where ipj = f{xj) and j = l,2,...,A^ — 1. On substituting the solution ipj = 
Q^jd Tf^Q obtain 

(27) 



Since (fij = e~^^^ is also a solution, then the general one will be (pj = Ae^^^ + 
Be~^^^ . From the boundary conditions y^o = = we obtain ipj = 2Ai sm{j6), 
where 9 = nn/N, n = 1,2, . . . , N — 1 (note that — )■ oo as /i — )■ 0). Therefore, 
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the approximate eigenvalues are 



1 — cos(2n7r/i) 
4/^2 



6" 



+ ... 



(2J 



This expression shows that the error in the numerical calculation of the eigen- 
values decreases quadratically with the step size h and increases with the 
quantum number n. A method is called k*^ order if its error term is of order 
^fe+i xhe naive shooting method just outlined is first order. There are 
other well-known numerical integration algorithms of greater order like the 
fourth-order Runge-Kuta method [12] that we will use in the calculations 
below. 

Tannous and Langlois [IJ do not indicate the value of Xq chosen for their 
calculation on the infinite square well. More precisely, they appear to be rather 
inconsistent about this important point. They first state that if the potential 
is symmetric in the interval [xi,X2] they choose xq = {xi + X2)/2. However, 
in the discussion of the boundary conditions, they say that if the potential is 
symmetric in an interval of length a they set a;i = and 0:2 = a/2 in which case 
we expect a different value of this coordinate point: < Xq < a/2. Therefore, 
in order to illustrate the application of the method to this trivial model we 
choose an arbitrary value < xq < 1 and do not exploit the symmetry of 
the equation about x = 1/2. The two linearly independent solutions that 
satisfy Eq. (jl]) are C{x) = cos[k{x — Xq)] and S{x) = sin[A;(x — Xq)], where 
k = \/2e. Obviously, in this case we do not have to bother about reaching 
constant values of l± (saturation) because the coordinate interval is finite. 
Upon substituting the Dirichlet boundary conditions into equations f jTBj) and 
f lTS]) (adapted to the finite interval) we obtain 

F{e) = /+(6) - L{e) = -k---^^^ (29) 

sm(A;a;o) sm[A;(l — Xo)\ 
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Fig. [U shows F(e) for < e < 10 and xq = 1/8,1/4,2/5. Note that F(e) 
vanishes at A; = nvr, n = 1,2,... thus giving the well known eigenvalues 
disregarding the chosen value of Xq. We also appreciate the effect of Xq on the 
form of the characteristic function F{e). 

Commonly, it is not difficult to derive approximate expressions for the con- 
vergent and divergent asymptotic forms of the wavefunction because they are 
straightforwardly determined by the asymptotic behaviour of the potential 
v{x). Therefore, it only remains to have sufficiently accurate expressions for 
C{x) and S{x) and their derivatives in order to obtain the eigenvalues by 
means of the equations developed in sections [2] and [31 This problem is easily 
solved by means of, for example, a suitable numerical integration method [12]. 
If y{x) stands for either C{x) or S{x) then such an approach gives us its val- 
ues at a set of points Xq — N^h, xq — Nih + h, . . . , xq, Xq + h, . . . , xq + N^h 
where Nl and Nr are the number of steps of size h to the left and right of 
xo, respectively. The number of steps should be sufficiently large so that y{x) 
reaches its asymptotic value at both xl = Xq — N^h and xr = Xq + N^h and 
h should be sufficiently small to provide a good representation of y{x). The 
numerical integration methods also yield the derivative of the function y'{x) 
at the same set of points which facilitates the calculation of the Wronskians. 
If the potential is parity invariant we only need to integrate the Schrodinger 
equation from xq = to X/? = N^h. 

In Sec. [3] we provided a rigorous basis for the SEM-CFM that is one of the 
goals of this paper. In what follows we illustrate those mathematical results 
by means of another exactly solvable problem. For concreteness we choose 

v{x) = (30) 

cosh (x) 

where —oo < x < oo. The allowed dimensionless energies are given by [81IT3] 
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e„ = -^(A - 1 - n)\ n = 0, 1, . . . < A - 1 

A=^(l + v/rT8^) (31) 

and the spectrum is continuous for all e > 0. It is clear that A — ?■ 1 as f o — !■ 
and there is only one bound state when 1 < A < 2 (0 < < 1). As fo increases 
more bound states appear. As a result there are critical values of the potential 
parameter for which e„ = that are given by the condition A„ = n + 1 or 
vo,n = A„(A„ - l)/2 = n{n + l)/2. 

Since lim|^|^oo ^(a^) = we have Rc{x) — )■ e"^^ and Rd{x) — )■ e^^, where 
k"^ = — 2e (we only consider the interval < x < oo because the potential 
is parity invariant). Consequently, the allowed energies are determined by the 
conditions 

W(Rc, C)+ = lim \C'(x) + kC(x)] e"*^^ = 

W{Rc, S)+ = jim [S'{x) + kS{x)] e-^"" = (32) 
for even and odd states, respectively. 

Since the potential fl30|) is parity invariant we integrate the Schrodinger equa- 
tion from Xo = to x/j = N^h. Fig. [2] shows the Wronskians W{Rci C) and 
W{Rc,S) and the ratios C{x)/S{x) and W{Rc,C)/W{Rc, S) for the arbi- 
trary values e = —1 and Vq = 2.5. We appreciate that the ratios C{x)/S{x) 
and W{Rc,C)/W{Rc, S) approach the same constant value as x — i- oo as 
proved in Sec. [21 Note that the latter reaches the limit at smaller coordi- 
nate values because the Wronskians take into account the asymptotic form 
of the solution explicitly. In other words, the WM requires less integration 
steps for the same accuracy. We also appreciate that x/? = 5 is large enough 
for the WM and SEM-CFM calculations in this case. In order to compare 
both approaches we find it reasonable to set h = 0.01 and N^i = 500 in the 
fourth-order Runge-Kutta method [12] built in the computer algebra system 



Derive (http://www.chartwellyorke. com/derive. html) that we use in all our 
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calculations. The numerical integration routine for the WM and SEM-CFM 
is identical and the only difference is given by the functions that we choose 
for the construction of F{e). 

Fig. 2 shows the Wronskians (!32l) and the ratios C (xr) / S{xr) and S{xji) / C (xr) 
for — 10 < e < and vq = 10. We see that both approaches yield the exact 
eigenvalues marked by black squares in the same figure. The main difference is 
that the Wronskians change much more smoothly than the canonical functions 
that cut the abscissae axis sharply at the eigenvalues. 

In many cases it is not difficult to derive the asymptotic form of the con- 
vergent and divergent contributions to the wavefunction. As another example 
consider the anharmonic oscillator v{x) = V2x'^ + f4X^ (f4 > 0). A particu- 
lar case is given by the double well discussed by Tannous and Langlois [Ijj. 
If we introduce the ansatz ^(x) = e~-^^^^ into the Schrodinger equation and 
keep only the dominant term we conclude that Rc{x) — )■ exp (— -^^-^3;^^ and 
Rd{x) — )■ exp (-^^-p^a^'^) as X — > oo. On inserting these asymptotic expansions 
into the WM or the SEM-CFM equations the saturation should appear at 
smaller values of the coordinate as illustrated above by means of the exactly 
solvable model potential f l30|) . Another interesting example is provided by the 
radial part of the dimensionless Schrodinger equation for a central-field po- 
tential v{r): 



where Z = 0, 1, ... is the angular-momentum quantum number and ^p{0) = 0. If 
we assume that limr^or'^v{r) = and insert the asymptotic behaviour ip{r) — )■ 
at origin we conclude that Lc{r) — )■ r'"*"^ and L^ir) — )■ are the regular 
and irregular contributions to the wavefunction. Taking into consideration 
these asymptotic behaviours in the WM or SEM-CFM equations will provide 
an advantage during the integration from tq towards origin. According to the 




/(/ + 1) 



+ t>(r) ip{r) = eip{r) 



(33) 



2r2 
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discussion above one expects that it will not be necessary to integrate too 
close to the origin to achieve saturation. On the other hand, the asymptotic 
behaviour for r — )• oo is determined by the form of v{r) as discussed above. 

5 Conclusions 

The main purpose of this paper is to show that the Wronskians provide a 
rigorous basis for the discussion of the SEM-CFM equations as well as con- 
cepts like saturation. We have seen in Fig. [2] that the Wronskians reach the 
asymptotic value or saturation at smaller values of |x| because they take the 
asymptotic behaviour of the wavefunction explicitly into account. However, 
the gain in numerical performance and efficiency derived from this result is 
not as important as giving the students the opportunity to discuss the asymp- 
totic behaviour of the wavefunction for a given quantum-mechanical problem. 
The most general equations for both approaches developed in sections [2] and 
|3] suggest that they are essentially identical. If, for example, instead of substi- 
tuting the boundary conditions in the SEM-CFM function F{e) we substitute 
the correct asymptotic behaviour of the wavefunction the SEM-CFM equa- 
tion for the eigenvalues becomes the WM one. In such a case that one finds 
it rather laborious to develop a suitable analytical expression for the asymp- 
totic behaviour of the wave function at some singular point (as, for example, 
in the case of the Lennard- Jones potential then it may be more conve- 
nient to resort to the raw boundary conditions satisfied by the wavefunction 
as discussed in Sec. |3] (for example, Eq. f l23|) ). In other words, we have the 
chance of using the SEM-CFM directly or improving it by explicitly using the 
asymptotic behaviour of the wavefunction. 

In closing we want do discuss some points that may be of pedagogical value. In 
the first place Tannous and Langlois [1] state that the characteristic equation 
F{e) = is not a matching condition but a dispersion relation. The matching 
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methods typically integrate the Schrodinger equation inwards from left and 
right and require that (yo)- = (yo)+ and (i/'q)- = {yo)+ at an intermediate 
point Xq. Obviously, satisfying these equations is equivalent to obtaining the 
roots of F{e) = /+(e) — /-(e). One can also integrate the logarithmic derivative 
y'{x)/y{x) inwards and obtain exactly the same equation. On the other hand, 
the SEM-CFM integrates the canonical functions (or normalized solutions) 
left and right from a given point Xq and then match the logarithmic deriva- 
tive for each case at that intermediate point. The difference between both 
strategies is merely the direction of the integration (inwards or outwards). In 
fact, Kobeissi [2] explicitly refers to the continuity condition of the eigenfunc- 
tion. The WM also proceeds outwards but does not focus on the logarithmic 
derivative. Instead it makes use of the Wronskians that become constant as 
the wavefunction approaches the asymptotic region to derive equations for the 
coefficients of the two linearly independent solutions. 

The authors mention that the SEM evaluates the eigenvalues directly and 
avoids losing accuracy associated with rapidly oscillating wavefunctions of 
highly excited states. One should not forget that the calculation of the canon- 
ical functions is equivalent to the calculation of the wavefunction. Kobeissi |2] 
explicitly shows the oscillatory behaviour of such functions for an excited state 
of the Morse oscillator. It is customary to state that the CFM does not calcu- 
late the eigenf unctions explicitly [THUE] when it is obvious that the approach 
already does it through the numerical integration of the canonical functions. 
In other words, the numerical integration of the Schrodinger equation explic- 
itly calculates two oscillatory functions but their ratios (or the appropriate 
Wronskians) do not reflect such oscillations as shown in Fig. [2j 

In closing we want to discuss the statement that the SEM enables 'full deter- 
mination of the spectrum in a single run' pLj. It is not clear what they mean 
because they calculate each root of F{e) = by means of the secant method. 
This algorithm requires the numerical evaluation of F{e) several times till e is 
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sufficiently close to the chosen root and each calculation of F{e) requires an 
outward integration of the normalized solutions from xq. In other words, just 
one eigenvalue requires many integrations. 

References 

[I] Tannous C and Langlois J 2011 Eur. J. Phys. 32 1519. 
[2] Kobeissi H 1982 J. Phys. B 693 693. 

[3] Tannous C and Langlois J, Marching toward the eigenvalues: The Canonical 
Function Method and the Schrodinger equation. larXiv:1003. 0184V 1 [quant-ph]. 

[4] Tannous C, Fakhreddine K, and Langlois J 2008 Phys. Rep. 467 173. 

[5] Leubner C 1983 Int. J. Quantum Chem. 24 127. 

[6] Kobeissi H 1986 Int. J. Quantum Chem. 29 163. 

[7] Tellinghuisen J 1988 Int. J. Quantum Chem. 34 401. 

[8] Fernandez F M 2011 Eur. J. Phys. 32 723 (|arXiv:1101.3209v l [quant-ph]). 
[9] Fernandez F M 2011 Am. J. Phys. 79 877 ()arXiv:1101. 09571 ^1 [quant-ph]). 
[10] Fernandez F M 2011 Eur. J. Phys. (in press) (|arXiv:l 107.409 2 v2 [quant-ph]). 

[II] Bonneau G, Faraut J, and Valent G 2001 Am. J. Phys. 69 322. 

[12] Press W H, Flannery B P, Teulosky S A, and Vetterhng W T 1986 
Numerical recipes. The art of scientific computing (Cambridge University Press, 
Cambridge). 

[13] Fliigge S 1999 Practical Quantum Mechanics (Springer- Ver lag, Berlin). 



16 




Fig. 1. Characteristic function F{e) for the particle in a box for xq = 1/8 (blue 
squares), xq = 1/4, (red solid squares) and xq = 2/5 (green circles) 
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Fig. 3. W{Rc,C){xr) (red solid line), W{Rc,S){xr) (blue dashed line), 
C{xfi)/ S{xfi) (red squares) and S{xr)/C{xr) (blue circles) for the potential ([30|) 
with vq = 10 
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